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Abstract 

We study the effect of coherent and incoherent broadening on neutrino os- 
cillations both in vacuum and in the presence of matter (the MSW effect). 
We show under very general assumptions that it is not possible to distinguish 
experimentally neutrinos produced in some region of space as wave packets 
from those produced in the same region of space as plane waves with the same 
energy distribution. 
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I. INTRODUCTION 



Neutrino oscillations have been the subject of intense theoretical and experimental re- 
search. To date there is no evidence for oscillations in terrestrial neutrino beams. The 
deficit of solar neutrinos can be explained by neutrino oscillations with a (mass) 2 difference 
Am 2 ~ 10~ 6 eV 2 together with the enhancement of these oscillations as the neutrinos pass 
through the sun by the MSW effect Q. There are also hints of possible neutrino oscillations 
with Am 2 ~ 10 _2 eV 2 in atmospheric neutrino experiments B. 

Approximately twenty-five years ago an interesting suggestion was made for probing even 
lower values of Am 2 using the 3% annual variation of the earth-sun distance ||. In this 
case, rather than simply observing a net average decrease in the electron neutrino intensity 
by an amount sin 2 (2#)/2 one could observe the actual oscillations of the electron neutrino 
flux. The idea is to use the z/ e 's from capture on Be 

7 Be + e~ -J Li + v e (1) 

which results in a neutrino energy E v ~ .86 MeV with a small energy spread. Thus if the 
neutrino oscillation length 

£osc = -7 2 (2) 

Am z 

is within one or two orders of magnitude of the variation AR ~ 5 x 10 11 cm of the earth-sun 
distance then, depending on the value of sin 2 (29), it may be possible to see the neutrino 
oscillations provided Am 2 is in the range 10~ 9 — 10~ n eV 2 . 

One of the most essential ingredients in making the above scenario work is that the 
spread in energy AE of the neutrino "beam" is not too wide. This is especially true in 
this case since R/ AR ^> 1. If AE is too large then by the time the neutrinos arrive at the 
earth the oscillation patterns for neutrinos of different energies get sufficiently out of phase 
to wipe out any potentially observable oscillations. This results simply in a decrease of the 
total v e intensity by an amount sin 2 (29) /2. A coherence length L max is usually defined as 
the distance at which a neutrino of energy E has undergone one oscillation more than a 
neutrino of energy E + AE. This coherence length is given by 

AnE 2 / E \ 

max ~ (Am 2 ) AE ~ osc \Ae) { ) 

and the total number of complete oscillations will be 

N = Lmax U) 

iv max t V / 

-^osc 

Thus when AE/E is larger than about 1/30 we can no longer observe the oscillations and 
a narrow energy range AE is therefore required. 

The argument above assumed that the energy spread of the neutrino beam is incoherent 
in origin in the sense that it is due to slightly different energies of various neutrinos. The 
main origin of this energy spread AE is that the continuum electrons which are captured by 
the Be have an energy spread AE e ~ kT which translates into a similar spread AE V ~ kT 
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of the emerging neutrino energies. Another slightly smaller contribution to AE U = AE 
originates from the different Doppler shifts due to thermal nuclear velocities (relative to the 
line of sight) - the analog of the well known Doppler broadening in atomic spectroscopy [|J. 

"Coherent broadening" - namely the quantum mechanical spread SE of a single neutrino 
can also lead to the loss of the oscillation pattern ||. The well known natural line width in 
atomic spectroscopy: 

6E ~ T ~ (Tde^r 1 (5) 

is an example of coherent broadening. The finite lifetime Td 0C ay of the level interrupts the 
classical emission of the wave-train and limits the size of the wave packet 5x to: 

6x = cr decay (6) 

with the momentum (5P = SE/c) and the configuration space (5x) widths being inversely 
related to each other via a Fourier transform of a Lorentzian to an exponential. 

Another example of coherent broadening is the collisional broadening (also known as 
the "pressure broadening") of the neutrino line. It stems from the interruption of coherent 
emission by collisions of the emitting atoms. The corresponding wave packet size is given by 
an analog of Eq. but with Td eC ay replaced by ^collision - the effective time interval between 
"relevant collisions" . This (nuclear) collisional broadening effect has been extensively studied 
as the major contributor to the loss of coherence in neutrino oscillation. There have been 
various estimates of the strength of the effect leading to estimates of the size of 5x = ct co iii S ion 

1,1,1, 0. 

A third contribution to the coherent broadening which we believe is likely to contribute 
even more to the energy spread 5E of the neutrino wave packet is the small size of the 
wave packets of the captured electrons. Since the K electron ionization energy in Berillium 
E wn = Z 2 R y = 16-R y ~ 220 eV is small in comparison with the thermal kT ~ keV energy, 
the capture in reaction (|T|) is primarily that of continuum electrons. An electron wave packet 
of size S e will traverse the (point like) nucleus in a time 

St=- (7) 

V e 

where v e is the velocity of the electron. Because the weak interaction underlying the capture 
process (0) is local, the time available for the v e emission is 5t and the size of the outgoing 
v e wave packet emitted with velocity c (c = 1 in our units) will be 

Su = - (8) 
The thermal kinetic energy of a typical electron is \m e v 2 e ~ \kT . Thus 



!3kT 

.08 (9) 



It remains only to estimate the appropriate wave packet size 8 e to be used in Eq. (0). 
The electrons suffer many random collisions in the hot core which tend to localize the wave 
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function and reduce the wave packet size. If the only information available is that the 
electrons are in thermal equilibrium then 5 e is expected to be of the order of the thermal 
wave length: 

2tt 2tt 

S e ; 10 

m e v e y/3m e kT 

which then leads to a neutrino wave packet size 

2tt 



8E U 



6 x l(r 8 cm (11) 



This 5 U is smaller (and the corresponding incoherent broadening is larger) than all previous 
estimates. 

The three mechanisms described above all lead to the conclusion that neutrinos are 
emitted in the sun as wave packets with a rather small size 5 V corresponding to "coherent 
broadening" of the neutrino line by an amount 5E ~ 2tt/5 u . This coherent broadening also 
leads to the loss of the oscillation pattern || after a coherence length L coh which is precisely 
equal to the coherence length L max derived in Eq. @. This result can be derived technically 
by decomposing the wave packet into plane waves of energy E with a probability distribution 

P{E) = |^ (E) | 2 (12) 

and repeating the discussion leading to Eqs. (|3|) and (f|). This leads to identical conclusions 
but with AE replaced by the energy spread 8E given by: 

(5E) 2 = J dE P(E)(e 2 - E 2 ) (13) 

There is, however, a simple intuitive explanation for how the oscillations are lost in 
terms of the wave packet of the neutrino in configuration space ||. Consider an electron 
neutrino wave packet which is emitted at t — from the solar core. At t — the v e can 
be written as a superposition of two wave packets with identical shape corresponding to the 
mass eigenstates \v\) and \v-i). 

\v e (t = 0)> = cos(0)h) + sin(0)|*/ 2 ) (14) 

This initial wave packet will quickly spread in the directions (x, y) perpendicular to the 
direction of motion but the spreading in the direction of motion (z) is negligible due to 
Lorentz contraction effects. Due to the different mass of the v\ and v 2 their wave packets 
travel with a different (group) velocity 

a Am 2 

/\ V = V2 - Vl = —— (15) 

Thus after a time t has elapsed and the neutrino has traveled a distance r ~ t from the 
source the two wave packets move with respect to each other by an amount 

Ar = Avt ~ — — r (16) 
2£/ z 
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Neutrino oscillations are simply the "beating" of the two wave packets as they slide relative 
to each other by Ar = A with 



A 



E 



(17) 



the wavelength of the neutrino. The oscillation length of Eq. (0) is then recovered as: 



L, 



osc 



{value of r for which Ar 



A} 



Av 



Am 2 



AttE 



(18) 



The total number of possible neutrino oscillations is simply the total number of wavelengths 
within the wave packet, A max = 5 U /X = E/5E. After this number of oscillations the two 
wave packets do not overlap at all and all oscillations are lost. Thus the coherence distance 
L co h which is the maximum distance over which we see oscillations is given by 



which is precisely the result of Eqs. (|3|) and (|j) for the case of incoherent energy broadening. 
Indeed once Ar is greater than the size 5 V of the wave packet the v\ and the z/ 2 will have 
completely separated spatially. We would thus expect that they will not interfere when 
interacting locally with an electron or nucleus in a detector. 

The main aim of this paper is to study whether the two effects discussed above namely 
the incoherent versus the coherent broadening can be distinguished. They are clearly distinct 
physical phenomena which can be controlled (at least in principle) at the source. In an 
Atomic Physics analog the Doppler broadening can be controlled relative to the natural 
line width by adjusting the temperature of the system or by confining the atoms to a 
narrow channel transverse to the line of sight [|J. The more interesting question is: Can we 
distinguish these effects at the detector? In this paper we shall show that in all physically 
interesting situations the answer is "no" . We shall discuss some simple cases in which this 
answer is clear and then we shall prove some general theorems which will show that under 
a wide variety of physically attainable situations these two effects cannot be distinguished. 



Our goal in this section is to see whether one can distinguish an incoherent ensemble of 
plane waves with a mean energy E and an energy spread AE from an ensemble of wave 
packets each with the same mean energy E and the same energy width 5E = AE. Before 
proceeding we should make one point clear. Even in the "incoherent" case in which we 
have an ensemble of plane waves these waves certainly do not have an infinite extent in the 
z direction (the direction of motion). In fact even if we took each "plane wave" (with an 
energy in the MeV range) to have an energy uncertainty of the order of 10~ 5 eV (which is 
certainly a great underestimate for the solar neutrino case) the corresponding wave packet 
would still be only of the order of a cm in size!! Thus when discussing "plane waves" we are 
in fact referring to wave packets which are much larger than those discussed in the case of 
coherent broadening but much smaller than any macroscopic scales in the problem. 




(19) 



II. COHERENT VERSUS INCOHERENT BROADENING 
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A. An Example 



Our aim will be to show that the two broadening effects discussed above cannot be 
distinguished. We begin with a concrete suggestion for distinguishing these effects and we 
then show what goes wrong with this suggestion. 

Let us suppose that we were able to measure the energy of a neutrino with a precision 
e which is much better than 5E = AE. We then expect that for an incoherent beam of 
neutrinos with energies in a range AE about E we could recover the oscillations by measuring 
the neutrino energy to the precision e <C AE. By plotting the observed neutrino count as a 
function of 



r' = r - 



(20) 



we should see oscillations up to a new distance 

L e = (E/e)L osc > L max (21) 

with no loss of statistics. Note that AE is replaced by e in Eq. (|3|). 

If, on the other hand, we began with a wave packet with energy spread 5E then, at a 
distance larger than L co h (Eq. (|19"D), the wave packet of the v\ and the z/ 2 are completely 
separated and one might naively expect that there will be no oscillations even if the energy 
could be measured more accurately. 

This argument turns out to be wrong and we can understand what goes wrong in a very 
intuitive way. If we choose to measure the energy very accurately (to an accuracy e) we 
require a time t ~ 1/e to make this measurement. If, during this time t, the second wave 
packet arrives at the detector then we will once more see the oscillations. The condition for 
recovering the oscillations is therefore 

Aw? 

t ~ 1/e > Ar ~ —r (22) 

where Ar is the distance between the wave packets and r is the distance from the source 
(Eq. flni|)). The oscillations thus persist up to a new distance L' e which is the value of r for 
which Eq. (^) breaks down. 

L' e = (E/e)L osc (23) 

which is precisely the same as the result ( |21| ) obtained for the incoherent neutrino beam. 

This behavior of the coherent beam is analogous to what occurs for a high Q oscillator 
hit by two successive pulses. The first pulse (in our case the v x beam) comes along and sets 
the oscillator in motion. It then continues to oscillate for a time t ~ 1/e during which time 
the second pulse (in our case the v 2 beam) arrives and causes the oscillator to be further 
excited. In this way coherence is maintained between the v\ and the z/ 2 beams even when 
they are spatially separated. What happens is that the accurate measurement of the energy 
picks out the plane wave in the wave packet which has existed coherently through both 
pulses. 

Our main goal will be to understand how general the above result is. In other words, 
to what extent is it true that an ensemble of plane waves will give the same result as wave 
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packets. Although there were some initial attempts to distinguish these processes it is now 
widely believed that they are indistinguishable. Our goal in this paper is to prove some 
theorems which clarify the conditions under which this is true and to show how general the 
result is. 

B. Measuring Observables which Commute with Momentum 

Before discussing the most general situation we review here the proof that the coher- 
ently and incoherently broadened neutrino beams lead to the same total rate and energy 
distributions for both i/ e 's and z/^'s. 

1. Oscillations in Vacuum 

Let us consider two cases representing two possible electron neutrino beamsQ leaving 
some region of the sun at time t — 0. In case a we have an incoherent mixture of neutri- 
nos each of which is a nearly ideal plane wave (with some extremely small energy spread 
SE pw <<< 5E). In this mixture the probability of finding a neutrino of energy E is given 
by some probability distribution P(E) which is centered about some energy E with a width 
AE. In case b all the neutrinos come with the same quantum state. This state is a wave 
packet with amplitude ^f(E) for a plane wave component of energy E. We choose this am- 
plitude so that the probability distribution |\I/(.E)| 2 precisely matches the distribution P(E) 
of case a . Consequently the widths of the two distributions are also equal: 5E = AE. In 
this section, for simplicity, we shall treat the plane waves of case a as ideal plane waves with 
6E pw = 0. 

At t — the wave function for the case b is given by: 

l^ = 0)> = £a p >,*> (24) 

p.i 

where the sum (which is actually an integral) is over momenta p in the z direction (the 
direction of motion) and over mass eigenstates % = 1, 2 and the a Pi i are chosen to give an 
electron neutrino with the appropriate wave function at t — 0. Since the \p, i) are eigenstates 
of the Hamiltonian, at a later time t, the wave function is given by 

\m) = E a P^~ le(;h \p^) ( 25 ) 

p,i 

where = \Jp 2 -\- mf is the energy of Vi with momentum p. 

Suppose now that at time t we measure an observable Q which commutes with the 
momentum operator. Q may, for example, be the total number of electron neutrinos in some 
range of momenta. This is, in fact, the most common kind of measurement which can be 



They could, of course, be any linear combination of electron and muon neutrinos. Electron 
neutrinos were chosen for definiteness only. 
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made. In this case Q has only diagonal matrix elements in momentum space. Therefore the 
expectation value of Q at time t is given by: 



W)\QMt)) = Y, 



(26) 



The expression inside the square brackets is precisely the expression for the expectation value 
(Q) p of Q for a plane wave which has a total weight (i.e. normalization) |a p ,i| 2 + |«p,2| 2 and 
a relative amplitude a Pj i and a Pt 2 for v\ and u 2 respectively at t = 0. Thus 

(Q) = E(Q)p (27) 

p 

which is precisely the result one obtains for the incoherent beam of case a . Thus the mea- 
surement of any observable which commutes with momentum yields the same result for 
case a and case b . 

Although this result may seem trivial it is in fact rather powerful. ^From this result we 
can verify the result claimed in Sec. [I] that if we use the variation in the earth-sun distance to 
look for oscillations in the neutrinos from 7 Be both the coherent and the incoherent neutrino 
beams give the same oscillation pattern. A-priori the above theorem is not applicable 
since the experiment involves measuring the spatial dependence of the neutrino flux which 
involves the use of an operator which does not commute with momentum. This is however 
an example for which the conversion of spatial to temporal dependence can be done reliably. 
Thus, although we measure the spatial variation in the neutrino flux, we can compute the 
temporal dependence of this flux by computing, for example, the total number of electron 
neutrinos with a given energy as a function of time. This estimate will be reliable since, as 
discussed previously, even the "plane wave" packet is still extremely small (certainly much 
less than a cm in size) relative to the relevant astronomical scales. 



2. Oscillations in Matter 

The above proof that an ensemble of plane waves cannot be distinguished, at the detector, 
from an ensemble of wave packets with the same energy distribution can be extended to the 
case of neutrino oscillations in matter (the MSW effect) ||. To this end imagine that at 
t = an electron neutrino is produced (at the origin) in matter in which the density of 
electrons (along the direction of motion of the neutrino) is given by p e (z). (This is of course 
an approximation in which we neglect variations of the density in the transverse directions.) 
The "vacuum eigenstates" and [1/2) are no longer eigenstates of this system. Instead one 
can find new eigenstates of the Hamiltonian which include the full spatial variation of the 
density. These eigenstates will of course no longer be momentum eigenstates. For relativistic 
neutrinos one should, in principle, solve the Dirac Equation but for the present discussion 
since spin is not a crucial variable it suffices to consider the Klein-Gordon equation |T0|| : 

(-(IH ,+ *)(£)-<;) (28) 



s 



where Mo is the vacuum mass matrix and the matrix A m accounts for the effect of charged- 
current scattering of the v e off the electrons in the medium: 

r , f m\ cos 2 9 + m\ sin 2 6 (m 2 — m 2 ,) sin#cos6>\ . ( \^2GFPe{z) 



M ° v (m 2 -m 2 )sinflcosfl m\ sin 2 + m 2 cos 2 ) ' Am { o" oj ^ 

The eigenstates of this system with energy E will no longer be eigenstates of p z (p x and 
p y are assumed to be zero) but will be labeled by some other parameter which we call 7. 
We shall call these eigenstates 

17,1) and | T , 2) (30) 

In the regions of space where the density vanishes these eigenstates will behave as plane 
wavesf] with some momentum p V ac(7)- They will correspond to vacuum mass eigenstates of 
the system. In a region of space in which the electron density is nonzero but nearly constant 
the eigenstates will again be nearly plane waves but now corresponding to the usual 
neutrino eigenstates in matter. 

Suppose now that at t = we prepare an electron neutrino in a state described by some 
rather narrow wave packet \ip(t = 0)). (This is case b of Sec. [II B 1| .) At t — 0, in analogy 



with Eq. (pip, this state can be expanded in the eigenstates \ j,i) described above 

lV>(* = 0)) = £«7,il7,*> (31) 



7,1 



We now allow the state to propagate to a later time t . At this later time the state is given 
by 



l^)>=E«7,e-^|7,*> (32) 



where eW is the energy of the state I7, i). In any reasonable case the size of the wave packet 
at time t will be much smaller than the scale of variations in the electron density. (This 
is especially true if the measurement is made in vacuum.) Thus Eq. ([32]) amounts to an 
expansion in the momentum eigenstates of the neutrinos in matter with density p equal to 
the density at the location of the wave packet. Every 7 corresponds to some momentum 
p(j) which depends on the density p. Thus 

toK*)) -2«p,<b>*;p) (33) 

p,i 

where for any given value of p and the corresponding value of 7 the coefficients d p i and d p 2 
are linear combinations of a 7i i and a 7j 2- 

Now suppose that at time t we measure some operator Q which commutes with (the 
z component of the) momentum. The off-diagonal matrix elements of Q vanish in the 
momentum basis. Thus Q will have only diagonal matrix elements between the various 



2 It may in fact be a superposition of an incoming and an outgoing plane wave if there is reflection. 
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\p, i; p) in Eq. ([33|). Since each of these \p, i; p) corresponds to one of the energy eigenstates 
it follows that the expectation value of Q is given by 



<tK*)IW(*)>~£ 



£<; e ' 



Mi 



(7>il Q 



— lei, t 



7,i) 



(34) 



This expression is ciiicilci gous to Eq. fl26|) in Sec. |11 B 1| . Each term in the sum is precisely 
the result which we would have obtained for the expectation value (Q) 7 of Q for a state 
which was initially in an approximate momentum eigenstate corresponding to 7 but with 
total weight |a 7i i| 2 + |a 7i 2| 2 and a relative amplitude a 7) j and a 7i 2 for v\ and z/ 2 respectively. 
(Recall that the realistic plane waves are actually extremely narrow on the scale of the 
density variations.) Thus 



(Q) = £<Q>- 



(35) 



which is precisely the result one obtains for the incoherent beam of case a . There is thus no 
difference between the wave packet ( case b ) and the plane wave ( case a ) ensemble even in 
matter when only operators which commute with momentum are measured. 



C. Unrealistic Measurements which CAN Identify Wave Packets 

^From the above proof it seems that a keen measurement which combines a measurement 
of both position and momentum information might be able to distinguish an ensemble of 
wave packets from an ensemble of plane waves. The simplest way of doing this would, 
however, require precise knowledge of the point of origin and the time of origin of the wave 
packet. Suppose, for example, that we knew that all the wave packets in our ensemble 
( case b ) were centered at the origin (z = 0) precisely at time t = 0. Suppose also that in the 
alternative scenario ( case a ) we also knew that each (nearly ideal) plane wave (which is still 
a wave packet but with a much larger spatial extent than that of case b ) in the ensemble 
was centered at the origin at t — 0. Under these assumptions about our previous knowledge 
and by a careful timing measurement at the earth to determine the duration of the neutrino 
pulse we could distinguish the two cases. (In fact in case b there may be two separated 
pulses.) This scenario is, of course, totally unrealistic and we shall see below that if we 
allow for an uncertainty in the location of the initial packets it again becomes impossible to 
distinguish the two cases by any measurement at the earth. 

There is another scenario under which it is clearly possible to distinguish the two cases. 
Suppose we have a detailed theory for the production mechanisms of the two cases which lead 
to some different observable at the source. Suppose, for example, that the position or mo- 
mentum distributions for the two cases are expected to differ. Then clearly such information 
can be used to decide which mechanism is producing the neutrinos (or, more realistically, 
which mechanism dominates). However in the case of level broadening we have no such 
information. Both the energy and the position distributions are expected to be roughly the 
same. The question which we are asking is: Assume we are given two "sources" of neutrinos 
(or production mechanisms) with the same position (z) and momentum (p) distributions. 
Is it possible to tell by measurements at the detector which of the two "sources" produced 
these neutrinos? 
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D. General Theorem 



This question can be set up more precisely as follows: Consider the following two modi- 
fications of the scenarios case a and case b discussed above: 

In case A we have a nearly ideal plane wave which is actually a wave packet of a fairly 
large size Az. (Recall that Az will typically be much less than a cm!). We imagine an 
ensemble of such "plane waves" each of which has a nearly precise momentum (in the z 
direction) centered about po with a spread Sp. Assume that each plane wave has exactly the 
same spatial location. (This is precisely the case a above.) 

In case B we have an ensemble of wave packets. Each wave packet has a spatial size 5z 
which is much smaller than Az and a corresponding momentum spread 5p = 1/dz which 
is precisely equal to the Sp of case A . Up to this point this looks exactly like case b above 
except we now allow each wave packet in our ensemble to be, at t — 0, at a different spatial 
location. We assume that the wave packets are produced in precisely the same region Az 
in which the neutrinos of case A are produced with precisely the same z distribution.^ The 
two cases are shown pictorially in Figure I. 

All the above information is given to the experimenter together with the additional 
information that the z and p distributions for both cases are equal at t = 0. The question 
is: With only this information can the experimenter distinguish with any experiment the 
cases A and B above? 

Intuitively one might guess that the answer is "yes" . There should be some way to tell 
if we are dealing with wave packets or with (almost) plane waves! But in fact the answer is 
"no"!! No experiment can distinguish the above two cases. 

The most general proof of this statement would proceed as follows: 
Step 0. Choose values for Sz = 1/Sp <C Az and for the mean momentum po which were 
defined above. 

Step 1. Begin with an arbitrary (smooth) but fixed expression for the wave function of 
the nearly ideal plane wave of case A . The only constraints on this wave function will be 
that it is centered (say) at the origin, that its spread in position is (a fairly large) Az with 
a correspondingly tiny spread in momentum about some momentum p. Then consider an 
ensemble of such states each with a different momentum p. Choose an arbitrary but fixed 
distribution for these momenta. The constraint on this distribution is that it is centered 
about the momentum po with the given width Sp. 

Step 2. Construct the Density Matrix for the ensemble described in Step 1 above. 

Step 3. One must now prove that it is always possible to construct the following, seem- 
ingly completely different ensemble, which, nonetheless, yields a density matrix identical to 
the one obtained in Step 2 above. We first construct a wave packet which is centered at some 
location z. We are free to choose the form of the wave function with the only constraint 



In a realistic situation both the "plane waves" of case a and the wave packets of case b will 
be distributed over a region of space much larger than Az. In both cases this excess spread is 
incoherent. It is thus sufficient to prove our result for the case when the wave packet is distributed 
in z by the size Az of the plane wave of case a . 
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that its spread in position is approximately equal to d~z <C Az with a corresponding momen- 
tum spread 5p = l/5z. We then construct an ensemble of such wave packets and choose a 
distribution of locations z with the only constraint that this distribution be centered at the 
origin with a spread in position approximately equal to Az. 

The claim is that we can always choose the distributions in Step 3 so that the density 
matrix for Step 3 is identical to that of Step 2. This then implies that any measurement at 
all which is done on the two ensembles at any time t gives the same result!! We also claim 
the converse of this theorem namely that given a "wave packet" ensemble constructed as in 
Step 3 it is always possible to find a "plane wave" ensemble as constructed in Step 1 with 
the same density matrix. 

Note how the mass eigenstates V\ and vi appear nowhere in the above discussion. The 
reason for this and, in our opinion, the power of this proof is that it relies entirely on 
properties of the system at t = at which time the state is a pure v e state. 



1. Illustration in the Simplest Case 

We can show the essence of the proof by the following simple example. We model the 
wave packet (of case B ) by a superposition of only two momentum eigenstates \p\) and \p2)- 
In Step 1 above we imagine having the state \pi) with probability \ot\ 2 and the state \p 2 ) 
with probability |/?| 2 ; (|a| 2 + \f3\ 2 = 1). The density matrix for this system is simply 

\a\ 2 \pi)(pi\ + \P\ 2 \ P 2)(P2\ (36) 

For implementing Step 3 we may construct an analogue of wave packets at two different 
locations as an ensemble consisting of these two states with equal probability: 

|^±>= a\ Pl )±P\p 2 ) (37) 

The density matrix in this case 

|^ + )<^+l + IV'-XV'-I (38) 

is precisely the same as the density matrix for case A in Eq. (|36|). 
This completes the proof in this simple case. 



2. Gaussian Distributions 

One case in which the Steps 0-3 above can be carried out explicitly is when all distribu- 
tions are Gaussian. Thus in Step 1 we choose the "plane wave" of momentum p to have a 
wave function 

\ p . p i ane) = 1 1 [ dl exp f-^f^) 10 (39) 

where a ~ 1/Az. We then consider an ensemble of these states with a Gaussian distribution 
of momenta p 
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(p-po)' A 



^\--^r) (40) 



where a p ~ dp = l/8z ^> a. The density matrix for this case ( case A ) is given by 

PA= (^ P ) jpsjJ"" dpew {-—} s 5 i^j |,><(| (41) 

We now proceed with Step 3 corresponding to case B . Consider a wave packet with mean 
momentum p centered at some location z: 

1 f^^uz^f ( 1 ~poY 



\z\ packet) = - j / dl e'^exp [ - v j \l) (42) 



27T<T p V P 



We shall soon see that the correct choice for a p is 



a 2 p = a 2 p + a 2 (43) 

which is approximately equal to a p = l/8z as required. We then consider an ensemble of 
these states with a Gaussian distribution of positions z centered at the origin of the form 

fdz exp (-HfL ) (44) 



The correct choice for spread in position, a z , will turn out to be: 



° z a 2 (a 2 + a? 



This a z is approximately equal to 1/a = Az as required. The density matrix for this 
situation ( case B ) is given by 



Pb = 7 — r / — r / dl dl' dz 



( 2z 2 (i- Po y (i'- Po y 

exp 1 



r2 



n- Aa 2 p 4a 2 



(46) 

With the choices we have made for a p and a z in Eqs. ( pffi ) and (f|5|) it turns out that 
the density matrices and ps are precisely equal. The calculation is straightforward 
and most easily done by computing the matrix elements {1\pa,b\1')- in order to compute the 
matrix elements of pa only the integral over p must be done. This is a Gaussian integral. 
For pb only the integral over z must be done. This is simply the Fourier Transform of a 
Gaussian. The result is the same for p^ and ps and is given by: 

(l\PAMl ) = ^=TT exp ( — ^ ^ ^df- (47) 



We thus establish, for the Gaussian case, that the two ensembles are identical 
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3. General Proof 



In the Gaussian case described above we did not use the fact that 8z <C Az. In the case 
of a more general shape for the "plane wave" and the wave packet we shall present a proof 
which does rely on this approximation. We conjecture that it is possible to slightly modify 
the theorem ^ so that it will be valid for general values of 8z and Az but we do not have a 
proof at this time. 

We begin again with Step 1 for which we choose a "plane wave" of momentum p to have 
a wave function 



\p; plane) = J dl U{1 - p)\l) 



(48) 



where the function f a (l — p) has a width a ~ 1/Az. We then consider an ensemble of these 
states with a distribution of momenta p given by some function g ap ip — Po) with a width 
a P ~ 3p = l/Sz ^> a. The density matrix for this case ( case A ) is given by 

PA = Jdl dl' dp [f*(l>-p)f <7 (l-p)g <7p (p-p )] (49) 

We now proceed with Step 3 corresponding to case B . Consider a wave packet with mean 
momentum po centered at some location z: 

(50) 



|z;packet) = J dl e 1 z a ap (l — Po)\l) 

which has approximately a width a p . We then consider an ensemble of these states with a 
distribution of positions z centered at the origin given by some function h az (z) with a width 
a z which is approximately equal to 1 jo = Az. The density matrix for this situation ( case B ) 
is given by 



Pb 



dl dl' dz a* (I' -p )a a M - Pq)K z (z) 



dl dl' [a* ap (l'-p )a ap (l-p )h a (l - I')] \l)(l' 



(51) 



where h a (l — I') is the Fourier transform of h az (z) which has a width approximately equal 
to a. 

The requirement that the two density matrices are equal is now simply stated as: 

a* ap (l' -p G )a ap {l -p G )K{l -V)= J dp f:(l'-p)U(l- P )g„ p (p-pa) (52) 

(It is now clear why the theorem, as stated, cannot be true in general. Given arbitrary 
smooth functions g a and f a with the restrictions described previously it is certainly not 



4 The modification we have in mind is to relax the unnecessary restriction that the shape of the 
wave packet is independent of z. It is reasonable to consider an ensemble of wave packets all of 
which have the same width but with slightly different shapes. The same could be done for the 
"nearly plane waves" . 
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possible, in general, to find functions a a and h a which satisfy Equation fl52|) since the 
integral in (521) will not always factorize in the required form.) The result is however valid 



when the width o v of g ap is much larger than the width a of f a . If a <C a p and if the function 
g Gp is sufficiently smooth 

m'- P )f*(l-p)g* p (p-Po) = 
flil' - p)U{l - p)ylg„,{l-Po)yJg ap {V-po) + O (53) 

Thus 

/ dp fZ(l' - p)U{l - p)g ap (p - po) ~ yJg, p (l-po)yJg, p (l'-Po) * J dq f*MUl - V + q) 

(54) 

Thus if we identify the function a a with the square root of g a and the function h az {z) with 
the square of the Fourier transform of f a {p) then the equality in Eq. (E21) is satisfied to 
order 5z/Az as required. 



E. Consequences 

Although the result proven above is not entirely general it is sufficient for all cases of 
practical interest. The reason for this is that we have actually proven three things. The 
result that measurements which commute with momentum could not distinguish coherent 
from incoherent broadening was completely general and did not depend on the shape of the 
wave packet nor on its width. Secondly the proof that for Gaussian wave packets the two 
effects could not be distinguished with any measurement was also general and it did not 
depend on the width of the Gaussians. Thirdly our extension of the proof to arbitrary wave 
packet shapes was possible in the limit 5z <C Az. A practical attempt to distinguish the two 
mechanisms of broadening would likely begin with a theoretical calculation which assumes 
Gaussian wave packets for simplicity. Furthermore it would likely compare the wave packets 
to actual plane waves for which Az — > oo. We have shown that any such attempt is doomed 
to failure. We conjecture that the result is more general so that for an arbitrary shape of 
wave packet it is possible to find an ensemble of nearly plane waves which mimic its behavior 
exactly. 

An interesting corollary to the result proven in the previous section is that one cannot 
tell, on an event by event basis, whether one has a wave packet or a "plane wave". The 
proof is as follows: Suppose it were possible, on an event by event basis, to distinguish a 
wave packet from a plane wave. It would then be trivially possible to distinguish the cases 
A and B above since in one case we are presented with a plane wave and in the other case 
with a wave packet. In fact in just one event we would know with which case we are dealing. 
But, as we saw in the previous section we cannot do this since the density matrices for the 
two cases are identical. It follows that no such determination can be made on an event by 
event basis. 



This result does not contradict the recent work of several authors 11 on the ability to 



measure the wave function of a single particle via a "protective measurement" . There are at 
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least two requirements for such a measurement to be possible. The first is that the system 
needs an energy gap so that successive (soft) measurements keep the particle in the same 
state. The second requirement is that it is known a-priori that the system is in an eigenstate 
of the Hamiltonian. Thus, for example, it is in principle possible to measure the ground 
state wave function of a typical atom even on a single atom but, if we do not know whether 
the atom is in an eigenstate of the Hamiltonian or in a superposition of eigenstates then this 
cannot be determined on a single atom. An argument very similar to that in the previous 
paragraph can be used to prove this result. In our case neither of these conditions are 
satisfied. We do not have a gap and we certainly are not in an eigenstate of the Hamiltonian 
when we are dealing with a wave packet and/or we start with a pure flavor state such as a 



The theorem presented in the previous section also provides a general tool for under- 
standing how the size of a Quantum Mechanical wave packet affects physical results in 
various circumstances. In fact the style of our proof which relies on the initial properties 
of the system rather than on the details of its time evolution is extremely useful. There 
have been several instances in which either careless approximations or faulty logic have lead 
to conclusions which disagree with our very general result. To illustrate this point imag- 
ine, instead of using our general proof, that we evolve each of the two ensembles to a later 
time t and then compared them. We must, by our theorem, get the same density matrix 
for each ensemble. But in doing this calculation we might make several approximations to 
simplify the calculation. We might, for example, neglect the longitudinal spreading of the 
wave packet. It turns out that even when this spreading is negligible compared to the size 
of the wave packet it has a significant effect on the final density matrices and we would 
find significant differences between the two ensembles. We know from our theorem that this 
cannot be the case. Indeed when the effect of longitudinal spreading is included all results 
computed with pa and p# agree. 



The main focus of this paper was the question of our ability experimentally (even in 
principle) to distinguish incoherent broadening of a neutrino line (such as the 7 Be solar 
neutrino line) from coherent broadening of such a line. Of particular interest was whether 
these two types of broadening would have different effects on neutrino oscillations and the 
MSW effect. We began by identifying processes which contribute to these mechanisms of 
broadening. Coherent broadening results from several processes including the natural width 
of the emitting nucleus, pressure broadening caused by collisions of this nucleus and the finite 
size of the wave packet of the captured electron. We argued that this last process leads to the 
smallest estimate for the spatial size of the neutrino wave packet (~ 6 x lCT 8 cm). Incoherent 
broadening results mainly from the thermal energy spread of the captured electron as well 
as from the Doppler shift due to the thermal motion of the emitting nucleus. 

We then began to present our argument that although the two forms of broadening were 
distinct physical processes which could be controlled at the source they could not be dis- 
tinguished at the detector. We first showed that if the detector had an excellent energy 
resolution not only could oscillations due to an incoherent ensemble of (nearly) monoener- 
getic neutrinos be restored but oscillations of a coherent neutrino beam could also be restored 




III. SUMMARY AND CONCLUSIONS 
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despite the physical separation of the v\ and the z/ 2 at the detector. We then proved that 
the measurement of any operator which commuted with momentum could never distinguish 
a wave packet from a plane wave. We extended the proof of this result to the case in which 
the neutrino propagates in matter (the MSW effect). 

The next stage was to show that if we had no a-priori knowledge of any difference in the 
properties of the coherent versus the incoherent neutrino "beams" there was no measurement 
which could distinguish them. Our method was to show that it was possible to construct two 
ensembles, one corresponding to "nearly plane waves" and the other to wave packets which 
had the same density matrix at t — 0. This would imply that the density matrices were equal 
at all later times and that no measurements could distinguish the two cases. We presented a 
complete proof in the case of Gaussian wave packets by showing that the density matrix at 
the source for an ensemble of plane waves with a given (Gaussian) energy distribution was 
equal to that of an ensemble of wave packets each with a much narrower z distribution but 
distributed, incoherently, over the same range of positions as the "incoherent" ensemble. We 
extended this proof to the case of non-Gaussian wave packets in the limit that the spatial 
size of the wave packet was much smaller than the spatial size of the "nearly plane wave" . 
We conjectured that the result is even more general and that given any ensemble of "nearly 
plane waves" with a given energy and position distribution we can construct an ensemble of 
wave packets which has precisely the same density matrix. 

There have been claims in the literature that wave packets could give different results than 
plane waves with the same momentum distribution. These differences show up either when 
the neutrinos are nearly nonrelativistic or when their momentum distribution is extremely 
broad so that 5p ~ p. This of course implies that some of the components of the neutrino 
wave function are nonrelativistic and that some of the neutrinos are moving "backwards". 
In all these cases it is essential to include the longitudinal spreading of the neutrino wave 
packet and to remember that if one calculates the number of neutrinos which should be 
observed at some location z one must compute the flux of neutrinos which involves the 
neutrino velocity. If these cautions are kept in mind one confirms the results of our theorem 
that there are no differences between the two scenarios. 

Although we have chosen to focus this paper on neutrinos and neutrino oscillations it is 
clear that the result is much more general. It applies to any particle for which the question of 
the distinguishability of a wave packet from plane waves is relevant. Some examples include 
neutral Kaon oscillations and the effect of wave packets in scattering theory. 
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FIG. 1. Pictorial representation of Case A and Case B described in Sec. IID 
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